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ABSTRACT 


The probability distribution of the minimum of the dis- 
tances from a randomly occurring trouble spot to n carriers 
patrolling along a line of length L is analyzed. Two approxi- 
mations, a Poisson process and a fixed lattice of equally 
spaced points, which bound the analytic model are developed 
and their usefulness and limitations are discussed. It is 
hypothesized that a two dimensional Poisson field and a two 
dimensional lattice of fixed points will form upper and lower 
bounds for the more algebraically tedious two-dimensional 
area patrol model. The hypothetical bounding distributions 
are developed for n units patrolling an area. Finally a 
one dimensional radius of influence model is developed which 
quantifies the contribution that the effective operational 
radius of a carrier airwing makes to the initial minimum 


distance analytic model. 
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ee neo OUeT LON 


The Commander in Chief Pacific has the basic responsi- 
bility of exercising control of the military forces of the 
United States from the west coast of North America to the 
east coasts of Asia and Africa. Since most of this vast 
eer Ener Ocean Or accessible by Sea lines of communi- 
@erron, the U.S. Navy has a critical role in the maintenance 
of a U.S. military presence in the Pacific Command area. 

The problem addressed in this thesis was posed by Vice 
Admiral E.P. Aurand, USN to Karl Eulenstein, Senior Analyst in 
Section J 77 at CINCPAC. As COMTHIRDFLT, Admiral Aurand 
wanted to know how many aircraft carriers would be required 
to patrol the western littoral of the Pacific Ocean Basen 
as well as the littorals of the Indian Ocean Basen in order 
that at least one patrolling carrier would be within a 
command specified distance of a randomly occurring trouble 
spot with at least a command specified probability. More 
formally, what is the probability distribution of the minimum 
of the distances between n carriers and a single random 
trouble spot somewhere along the patrolled coast 


This thesis addresses the analytical expression of, 


Pld = min(D,,D5,---,D,) <d] (1) 


10 





where: 


n is the number of carriers on station 


11P2 Pees Dd. are the distances from each carrier to 


a random trouble spot 
d is a command specified distance measured from 
Ene tCheeupLe SDOE- 
The model allows the determination of the minimum n required 
in expression (1) such that the probability 1s at least as 
large as p. Other related inferences and applications are 
discussed. 

In the following a one-dimensional model and two approxi- 
mations that give useful bounds to it are developed. Then 
the analogous two-dimensional bounding models are considered. 
Finally the concept of a "Sphere of influence" is incorporated 


into the one-dimensional model. 


Ile 





II. INITIAL ANALYTIC SOLUTION 


A. ASSUMPTIONS 

To increase the tractability of the analysis several 
Simplifying assumptions were made: 

(1) The complex curvilinear coastlines which mark the 
western boundry of CINCPAC'S area of responsibility are 
approximated by a Straightline of length L. 

(2) Each of the n carriers were asSigned individual 
patrol line segments of length = miles. 

(3) A carrier's position is assumed to be uniformly 
distributed on its asSigned patrol line segment. 

(4) A trouble spot will appear at random along the 
coast line. 

(5) An individual carrier is considered as a point. 
(The effects the carrier airwing's radius of influence are 
analysed in Chapter V.) 

(oy Each carrier's location within its patrol segment 
1s independent of the locations of the other n-l carriers. 
The location of a trouble spot T is independent of the loca- 


tions of the carriers. See Figure l. 


Eee DEVELOPMENT 

The Mean Value Theorem for Integrals is utilized in the 
derivation of the analytical results detailed below. fThis 
theorem shows one may approximate a probability density 


PametLOn at a point h as 


2 





Schematic and Symbology for the Analytic Model 


12 
A 
) K—  D, =. & 

— Da 
Dy + 
Ly —" | 
X, Xa X3 pe Xn-i Xn 


bem. = 1,2, «+e, N ERewloOcarcion OL Carrier 1 within 
line segment 1 


the distance to carrier 1 measured 
from a randomly occurring trouble 
Ssjolenes 4s 


O 
a 
if 
= 
NO 
S 


L length of line to be patrolled 


n the number of carriers assigned to 
Patwoewear line of length L 


1 BoanGdemLy G6CeuLring trouble spot 


tNotationally this paper will follow the custom of using 
uppercase letters for random variables and lowercase letters 
for parameters and outcomes with the exception of the 
parameter L which will denote the length of the line to 
be patrolled. 


Figure l 


is 





f(h)Ah = lim P{X¢ (h,htAh) ] 


Examination of the geometry of the problem suggests 


that the development of the cumulative distribution function, 


e(a)~ = PD = min(D,,---,D.) < dl be segmented as follows: 


(1) Left end segment such that Te (QO, = ) 


(2) Middle segments such that te (Ey for some 


eee eos 6 sey oe 
(3) Right end segment such that re (21) 
| a ae Ma X m= ae : X 
O — (m-IL L 
n a 


Figure 2. Line segmentation 


The general analytic result was developed by conditioning 
on the segment containing the target. As a result of the 
assumption that the trouble spot will appear at random along 


the line (0,L), the target is located in segment i with 


tthe cumulative distribution function which will be 
developed is actually a function of the minimum of the 
distances, d; the length of the patrolled line, L; and the 
mmeer OF pDatrolling units, n;Fpj(d;L,n). For simplicity 
this terminology will be abbreviated as FY (d) . 


14 





eeebability l/n; that is, 


(i-1)n in L 
Pa Pecan 7 ae a ae met ob — 1,2; ree 
ez) 
Given that re (2,2), we further condition on T being 


located in the first or second half of this line segment. 


bei 2D. and 


i-l 7 


D = [min(D,_,, 


whose cdf is being developed. The conditional pdf of D 


D, = |T-X;|, D aie [oo 


tap i+1! 


ea) oe where D is the random variable 

given T 1S unconditional over T (being located in a specific 
half of a line segment) and integrated over the range (0,d) 
to determine the conditioned cumulative distribution function 


edt) of D = [min (D,_,,D;,D )}] given that T is an element 


i+l 
of a specific line segment i. 
1. Left-hand Segment 
The assumptions made in developing this portion of 
the model are: 
(lL) @~Uniform (0,L), so Te (0,2) with probability 
aig) 


L 
(2) X, ~ Uniform (ne) 





eeer ung -orm (=, —) 
fa 
(4) Xi Xo T are independent random variables. 


a. Conditional Probability Density Function 


Let S({i) be the event that the random variable 


Gea 
on 


a ) of the line; then 


ieoccurs in the ee segment ( 


fot) } = l/n. 


1S, 





Consider the case te (0,54), so 
DS DS eer 


and the maximum value of D is Ht, PoweOe< doe, 


faiz,s(i) (alt) 4a = Se 
= P(X, ¢ (t-d,t-d-Ad) or xy e (ttd,t+d+Ad) } 
= ele 
L 
Port <d< = -t, 
P[De (d,d+Ad)] = P[X, « (ttd,ttd+Ad) ] = <Adn/L. 


; lope 1S me ide 
Consider next the case te (sn-5). so D = min (D,,D.) and 


Peeee- For 0<d< =H, 


P[D ¢ (d,dtAd)] = P[X. ¢€ (t-d-Ad,t-d) 


il 


Or Xy eetetadr,era+rhd))|] = 2Adn/L 
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Pip ¢ (d,dtAd)] = P({X, e (t-d-Ad,t-d) and X. > t+d+Ad} 


2 


or {Xy < t-d-Ad and X (t+d,t+d+Ad) }] 


2 = 


_ YC oe n _q).my).(4dn 
= (Adn/L) ag e d)-Ft ((t d) 7) ( L ) 
= sa(B)*[SB-t-d-ad+ (t-d-ad) ] 
n 
Sova. fl = micros) 
Thus it 1s seen that: 
2 ie ORs. 
= IE 
Eni p (alt) = Por 2eae (0,5) 
n/L ; tsd< ()-t 
DG) L 
Ties ls 
= £Oxr te (Ga) 
Asjo| nd ee 
= Cl a te 


ay 





b. Unconditioning Over te (0,5) And 


ee ivenns (1) 


7a, Tyae' 


The conditional probability density function is 


graphed below, where values of f (d|t) are shown for various 


Das 


regions of the d-t plane. 


— 
a) 





O 
OQ 
| 
| 
} 
as 
nN 
Baoure 3. Unconditional probability density function 
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We now find the conditional distribution of 


given S(1), by taking the expectation of f ) (d]} t) 


io) (20,5 Cl 
with respect to the conditional distribution of T given 
S(1). Thus unconditioning over T (given S(1)) has the 


Mathematical form: 


L/n 
foisczy 9) = J foie, scay CAIt) fqyg cay (t) at 
0 
moe0) < d < L/2n, 
d = 
= n, 8 an. a 
mm) J Up ap Ee +t I i t.-* 
0 d 
ib 
- a d 
n n n 
+ ; if 7: a> z dt 
=~ — 
= Te ey oe n, 2 mene = Ee 
eee eo Te2t) (= 2d) + 2(-) eS 
Migezer 215 2 
= (>) a Ge ond, | 
For = eG) _ P 
n-~—n 
ula L 
> ae 9) ae Nee 
Eoys(1) = f Aenea = 
0 d 
= (By2 LL ny2(, _Ad) (LL 
a eet Gd) +2 (=) (1 .? = ac) 
= (2)? (22 - 5a + 24°n/L). 


IEG, 


D, 


Le) 





Therefore in the lefthand segment, 


(B)*(S2-a-2na*/ul; for o<a< = 
Sol) S (1) sey 

(2) * (34 - 5a 4 24°n/Ll ; for + 2 T6l< = 

(4) 


tacmeonGtelonal Cumulative DLStribution Function 


ane 
OL 


The -cumulative aistribution function (cdé£) 
D conditioned by S(l1) in accordance with expression (2) 1s 


calculated by integrating the pdf over the range (0,d) of 


ys 
qd 
“pies ee J fy) g(q) (mn) dh 
Hom 0 < 9d < ee ’ 
= een 
qd 
Ghee 2 
ois (1) @) Ge lees = 20/7) dh 
0 
Mines > ce d 2nd 
= Ga > a7] » and 
L meee 
Pp| aan? 24 © 


20 














d 
= mie ok n 19 
Fp oie) ae | Gy Sh + 2h- T] dh + 57 
iy 2n 
- (2)? (3bd _ SA", and” _ (3%) L Hy + (3) (2) * 
eS ie n 2 
Dig ee be, 
Sm On om 
= (Ry 2/3Ld _ 5d* , 2na* _ 2317 | 19, 
L n 2 3L Th a Pgs! 
24n 
For the left-hand segment, 
2 3 
neoeerdeide a. 2nd IL 
nen bee = 30 es OTSA C Sor 
e ee) 2 )3td 25a 2nd 23t*, g b e a c k 
Bs (1) 24° ‘L n 2 3L Dine en 
24n 
a nes Lorvd<- = : 
sacoaiic) @ ' 
(5) 


Zee Middle Line Segments 


We next develop the cumulative distribution functions 
for the minimum of the distances between a random trouble spot 
miemm=2 Carrlers assigned to patrol stations 1 = 2,3, ..., n-l; 
stations which have at least one other patrol unit on either 


Side of their patrol segment. 


Cx AL 





The middle line segment conditional cdf is identical 
for each of the "middle" patrol stations and is distinguished 
from that Ree incsa with either of the end patrol stations 
as follows: 


(1) end segments 


D=min(D,,D,) or D= min(D._,,D,)- 


(2) middle line segments 


D = min(D;_,,/D;,D;,,) EOm—~peiea = 2a 3.) 44a g NOL 
; ; , es 
A typical middle line segment 1S segment 2, ae) = 


The assumptions for developing this portion of the model 


are: 
eae, . are 
Ojai Wntnornn.(0,i), so T Sn a) wlth DProbabllley 
L 
oud 
OQ Se Sfpigoan tS 
a ao aes 
Ee 28 
(3) X5 Uniform Sas a 
ls 
(4) X. ~ Uniform Sel, 


(5) X,1X5/X3, and T are independent random variables. 


aeeeconadlelonal= -robabililty Density Function 


Consider a given te (2,34), SO 


oO 
il 


Jt - X, | ’ 


il 
D = min{D,,D,} ; 


Dee, 





: 2L iG 
and the maximum value of Dis mer = eo} aan) el ae meee 


ED) p,s(2) (alt) 4d meee d,atAd) | 

= P[X, e (t-d,t-d-Ad) or Xo e (ttd,t+td+Ad)] = 2Adn/L. 
Ron: 7 = —de< ae ’ 

mb e— (d,dtAd)}|] = P[{X, « (t-d,t-d-Ad) and X, > (ttd+Ad) } 


or {X. ge (ttd,t+td+Ad) and Xy < (t-d-Ad) }] 


_ dn, 2b _ ny 440 yp g-aay- 
ee i (t+d+Ad) ] 7; t+ =~ [t-d-Ad] 
2Adn* L 
= a oa = Ce ely 
L 
DO 2asclal n 
= ae (d+Ad) =] 
Consider nexta given te Enso D, = Pee. = | Eo XS 
Din erie 2 AS aaa oo 


= min(D.,D,) and the maximum value of D is t -2. 
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Bor Us- ad < 


P[D ¢€(d,dtAd) ] 


Z 


P[X, « (ttd,ttdt+Ad) Or X, € (t-d,t-d-Ad) ] 


2Adn 
— 


Ber 4h =e Gl, oe 
n — - n 


Eibe (d,dtAd)j] = 


Thus it is seen that: 


P[{X, € (t-d,t-d-Ad) and A. 


{X, e (ttd,t+d+Ad) and X. < (t-d -Ad) }] 


2 


Adn 3L Mes an Lin 
er mam (t+d+Ad) JF + = ft a=Ad oe 
Dive L 
aan d - Adj 
L 
26eR 1 (d+Ad)n/L] 
2n i, 
TT Ome Ces Ca 
1 
Zi dn L Pie 
pit = di anos a 
au ca 
Tr Ome < = te 
ic 
an nd i. ae 75 
tia. go aaa a eos 
(6) 
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beeUMcenaterening Over te (=, se) And 


te (32,2) Given S(2) 


The conditional probability density function is 


graphed below, where values of f (cue) are given for 


B) | 2 


various regions of the d-t plane. 


Nn 





meaure 4. Unconditional probability density function 


PicmeeicittlenalieatsSerlbucion Of D, given S(2), 


1s now found by taking the expectation of f (dey ween 


Dies (2) 


Ze 





respect to the conditional distribution of T given S(2). 
mits, Unconditioning over T (given S(2)) has the following 


Mathematical form: 


Daley 
fois (2) @) i Eo r,s(2) (alt) Emig (2) (t) ace 
L 
n 
Por 0 < ad < a 
= — 2n 
BES (db/aa ot _g 
2n gle = Ave Ze 
moys(2) ‘° J aoe? OE 
i/o d+— 
aM ieh 
2n grouene jg! 
+ f Tat -=) (7) dt 
(2L/n) -d 
eS Tale 2 nd Ie Mace L 
= 2 (+) Coe ee =) a 2(F) a eae) 
‘ane gic call Za 
sp 2(F) (1 - =) Co an" cl} 
Z 
= 2(2)*(2 - SBE), 
Kaye = < qo =< ay 
n— —n 
Foe, e 
n het, 7 al L n 
Ep} 5 (2) (@ i = lh a) lela _S Sa bated (Gs) che 
16 d+ (=) 
n 
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“oS (2) ee 2B) (1 - BS) (Ea - 
+ 2(2)%(1 - BB - a - 4) 
ghee aay Osi 2nd? 
= 2 (=) (aa CISt a 7 ). 
Therefore in any middle line segment, 

Meee: 2 2nd? L 

2p) (E - DA ies, en 

eis (2) °°’ = 
2(B)* (<2 4a 2nd, mstsg-: 7 


c. The Conditional Cumulative Distribution PURCEION 
The cdf of D conditioned by S(2) is calculated 


by integrating the pdf over the range (0,d) of D; 








d 
Fois(2) (@) f Eois(2) (A) dh 
For 0<d< <, 
— — 2n 
d 

eau 2m 

Fp js (2) (h) eee) (= =) dh 
0 
2(8)2 (Ea _ 2nd”, 
7 if n 3L 
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and 


Fpig (2) (@ 


Pog (2) @) 








5 
DS(2) 2 OO 
LG L 
on = © Se 
d 
5 Tie aoe Pq ion 
Aon 2) (on ese 7 
iL 
Zi 
5 Meals La 2 2nd 
e+ 2 (Fr) 7 2d + 3 
Pal ie ees 
2) rn 2 (By 4 2bd _ 542 i 2nd” 
6 L n cye 
the middle line segment, 
3 
n> eae 
1G n Si 
5 Meee ie Zee 2c 
Ae 


oe Right-hand Segment 


Development of the conditional cumulative distribution 


ea 
ZL iG dire, 2 
(ea ore pease 
IL? 
= 1 
en. 
fOG O<d<—4 
2 
Tey); for ge<d<= 
12n 
L 
Le peat e ee = 
(8) 


function. of the minimum distance for the right-hand segment is 


28 





identical in methodology and form to the results obtained 


in analyzing the left-hand segment in paragraph II.B.1 above. 


We therefore omit the details. 


fe UNCONDITIONAL CDF FOR n SHIPS ON A LINE OF LENGTH L 

The cumulative distribution functions derived thus far 
Memeo etonea son S(i), i = 1,2, ..+, n. Since T~ Uniform 
(0,L) and each patrol segment is identically = lp) Lies aepelor 


the probability of T occurring in any specific segment is, 


P{S.] = ptr, (22, L241), = p(fatbDb_ at 


te 
aL n n n A I/E * wana 5 


Therefore the unconditioned cdf is derived from (5) and (8), 


as follows: 
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minimum of the distances from T to n patrol units is: 
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III. APPROXIMATIONS TO ANALYTICAL RESULTS 


A. RATIONALE FOR APPROXIMATIONS 

The existance of a simple, readily applicable approxima- 
tion, especially one with well understood and widely known 
characteristics of behavior, might be useful to the military 
commander and planner. This thesis investigates two such 
approximations and evaluates the range of values for which 


each may be an appropriate approximation of (12). 


fee FOLSSON PROCESS 

The Poisson process as a mechanism for positioning n 
units on a line of length L is appealing because of the well- 
developed and tractable body of theory regarding the exponen- 
tial distribution of the distances between units. The limi- 
tations of the approximation will be discussed after the 
model has been completely developed. 

1. General Properties 

For the model let n be the average number of ships 

located on a line of length L, so A = = is the average den- 
Sity of ships per unit length of L. Let Ws be the distance 
Meceween unit 1 and i+1l, for all i: i=1,2, ..., n-l. Then 


for a general Poisson process one has: 


(xL)2 q At 
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The interarrival "distances" (W) between adjacent ships are 
independent exponentially distributed random variables such 


that 


Ws 8 lapel. ee 
In what follows only those cases where T occurs between 
two units will be considered (the interval (0,X,) will be 
discussed later). Let the trouble spot T occur at random 


on the line of length L as illustrated below. 


c—— W ——_4 
Xi t eed 

a) ft tA 
K— MI —K— M2 —5 


Bene, 


Let Ml represent the distance from x, to t, M2 the distance 

imom t to Keays then it is easily seen that Ml1+M2 = W 
meal in aco ; 

the distance from the 1 Wneee coo Ene (ls 1 UNLeE. ne 

memoryless property of the exponential distribution implies 


that both Ml and M2 are distributed as exponential (1/)) 


random variables. Now with D = min(M1,M2), we have, 


nel) 


P[Min(M1,M2) <d] = P(D<4d) 


1-P(D>d) = 1 - P(Ml>4d,M2 > 4) 


1 - P(M1 > d)P(M2 > d) 
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and since Ml and M2 are identically distributed, 


Fy(d) = 1- (1 - F,(a)1%. 
Thus one has, 
CO i SG Wah 
ci ae (13) 


The model developed above for approximating the 
distribution of the minimum distance from a’random trouble 


Spot to the adjacent patrol units implies the cdf. 


> 2nd 
So ic ee Oa aldo sSaerel (14) 
mien 1S an exponential 2A cdf. 
2. Limitations To The Poisson Approximation 


The Poisson approximation cannot guarantee that n 
units will be positioned within the (0,L) limit of the line 
to be patrolled. For example there exists a positive proba- 
bility, however small, that no units will occur within 
mo) . 

The model as developed does not consider the specific 
Seseributions involved when T occurs to the left of the 


location of the first unit. The effect of this simplification 


38 





is to increase the size of d for an arbitrary probability 
p which results in a weakening of the boundary condition 
developed in III.D. 

The Poisson approximation does not preclude three 


or more units being very close to each other. 


C. FIXED STATION MODEL 

A second approximation, the fixed station model, is now 
considered. This model assumes each patrol unit maintains 
a fixed station located at the midpoint of the assigned 


patrol line segment x = Zate IL by/ 219) - 


xX a X3 Xu 
SS SS ee 
° = it oe - 

‘+ 2 f 


fe ode! Develooument 


Since the distance between adjacent units 1S constant 
(=) and the distance from any unit to its patrol line segment 
boundaries is L/2n, the model developed for an arbitrary 
patrol segment will be representative of all segments. 

Consider segment 1. Given S(1), T is distributed 
Maacormly on (0,2) and the model being developed will specify 


meme COnditional cdf of D = min{=-T,T- 5+}. That is, 
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- _ L L 
Pois(1) Seco nee = Pie d|T <=) P(T< =) 
+ P(D< @/T>S)P(T> 4) 





= ae Lyi es yi 

ed ee SE = 5 S/T > 5-) 5 

- Thy es L iy 5 de 

oe Ie 

d+=~--=- 

me Jey L oe 2 2Th 

=) ote ely tte a7) 1 ot 5 Ton aS ] 

2 

ey ct 2n IL 2596) 

fll eet oy ara 

ee Ue Poe ae) nd 

eee ye 

Be ee eed od 

— 2 215 L 

eo Ags L 

= ane. Or 9s! SS oy CLS) 


PeoEesstom (tos) 1s sconditieoned by the occurrence of 
Sl). However the conditional cdf given S(i) for each 1; 
meme 2, ..-, Nn? Of D = min(D,,D;,,) Uouueenmercal to (15) 
and the summation of n conditional cdf's multiplied by the 


Beneitioning probability l/n results in: 


Poa L 
PD ad) = ee eS Sl KGalo ECS 5 (16) 
which has a median value of ot (cz) 





2. Limitations Of The Fixed Station Model 
Ti .ol oust tatlonwor the model 1S the loss of 
realism in constraining a mobile unit to a small geographically 
fixed station. AS a consequence of this constraint the 


) is 27 Pout he wElxea Station 


maximum value of D = min(D.,D. 
r+ 2n 


di 


model. 


D. COMPARISON OF ANALYTIC AND APPROXIMATE MODELS 

The easiest method of contrasting the two approximations 
with the original analytic model is through a tabulation of 
meiere individual cumulative distribution functions for 
Specific values of n, L, and D. In Table 1d is tabulated 
aS a percentage of the ratio =, This procedure allows one 
to readily compare the various models at differing L and n 
values. 

For the L = 30000, n = 6 case the distribution function 
for the analytic model is bounded below by that of the Poisson 
process and above by that of the fixed station model. 

Migire (5) 01S a graphic presentation of the "fit" of 
the 2 bounding models when compared to the analytic model 
with parameters L = 30000 and n= 6. The graph depicts the 
residuals obtained from differencing the analytic and 
approximation models at various values of dad. Data points 
are from Table l. 

For the case under consideration the Poisson model median 
value of D (0.3465) does not approximate the analytic model 


median (0.266) as closely as does that of the fixed station 
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model (0.25). Except for values in the neighborhood of 
d= .5 the Poisson process model residuals are of greater 
Magnitude than Bway Ene fixed station residuals. 

The Poisson model gives reasonable (if somewhat conserva- 
tive) estimates of the analytic model results within the 
Benge 0.0 <d< .3. 

For this case, the fixed station approximation model 
not only bounds the analytic model, from above, it also 
appears to be a reasonable approximation of the analytic 
model over most of the range of de (0.0,1.0). The fixed 
Station median value of 0.25 is very close to the analytic 
model median value of 0.266. Using a criterion that a 
residual value greater than 0.1 indicates that the approxima- 
tion model is inappropriate for that value of d, it is seen 
that the fixed station model is appropriate for all d except 
0.4 < d < 0.6. This result, while developed for the case 
n= 6 and L = 30,000 will hold in general except in those 
unusual cases where n > L. 

The results tabulated in Table l suggest that there 
exists a natural hierarchy or ordering of the three models. 
meme) fixes the P(D<d) = .5 and then ranks the models 
according to the relative magnitudes of d required by each 


model to obtain this probability the result is as follows: 


Model Median 
fol s Sol 0.3465 
22 eanalvtic Oo PASS 
3. Fixed station OAs 
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Arranged in this manner the models reflect the following 


ecosts’ : 
1. Increasing level of control by higher operational 
authority 
2. Stochastically decreasing distances to random 
targets 


3. Decreasing randomness of own units motion 

4. In a hostile environment an increasing risk of an 
adversary accurately determining the positions of 
the patrol units 

5. Decreasing probability of adjacent units being in 

close proximity (overlapping of coverage) which is 
paralleled by decreasing problems with mutual inter- 
ference and contact identification. 

The geometry of the residuals plotted in Figure 5 indi- 
cates that the average of the two boundary cdf's might form 
a more reasonable estimate of the analytic model than either 
of the boundary models does when considered individually. 

The boundary model data in Table 1 was averaged, sub- 
tracted from the analytic (L = 30000, n = 6) model and the 
residuals plotted on Figure 5 as the average model residuals. 
The resulting residuals are obviously a very reasonable 
conservative approximation and lower bound for the analytic 


megel, particularly for those de (0,.65). 
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IV. APPROXIMATIONS FOR TWO-DIMENSIONAL PATROL AREAS 


The previous chapters have developed a one-dimensional 
patrol model and two approximations which form upper and 
lower bounds for the analytic model. It seems a reasonable 
supposition that when the problem scenario is expanded to 
include units patrolling a two-dimensional area that a two- 
dimensional Poisson field would form a lower bound of the 
two-dimensional analytic model and that a two-dimensional 
lattice of fixed stations could reasonably be postulated as 
forming an upper bound on the two-dimensional analytic model. 

The algebra involved in developing the two-dimensional 
analytic model for n randomly patrolling units is quite 
involved and time constraints precluded its development in 


this thesis. However the two-dimensional approximations are 


tractable and are presented here. 


A. TWO-DIMENSIONAL POISSON FIELD 
1. General Theory For A Poisson Field 
Consider an array of points distributed over a 
plane P. For each subset B of P let N(B) be the number of 
points within B. The array is said to be distributed in 
accordance with a Poisson process with density ¢ if the 
following assumptions are met: 


(a) the number of points in non-overlapping areas 
are independent random variables 


(b) for any subset B of finite area, N(B) is Poisson 


distributed with mean tb, where b is the area 
Gt eBee|Ret. 1] 
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points 


The dimensional analysis of ¢ is Tength. squared: 


Thus if N(b) represents the number of points in a region 
Partch area b: 
Pts) ) = 1) = ’ (20) 


The probability that there are no points within a radius r 


of a given location is determined by 


P(R>r) = P(no points within a region of area nr’) 


Tare : (21) 


2. Poisson Field Model 
a. Assumptions and Symbology 
The following assumptions are made: 
(lime ine MeUnits are distributed stochastically 
by a Poisson seasane with rate C. 
(2) The trouble spot T 1s equally likely to 
occur anywhere within the patrolled area. 
(3) The area to be patrolled is L units in 
tength and W units in width (L>W). 
b. Development 
In accordance with the theory for a Poisson 
n 


field, ¢ = iw and following expression (21), the Poisson 


field model is, 
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c. Limitations to the Poisson Field Model 
The Poisson field model cannot guarantee that 
n units will be positioned within the patrol area. Nor can 
the model preclude five or more emalige from being assigned 


positions which are in very close proximity to each other. 


ims LXED LATTICE OF POINTS 
The cdf of the model which formed the upper bound of the 
one-dimensional analytic model was a series of patrol 
Stations fixed in the middle of each patrol line segment. 
A reasonable two-dimensional approximation would seem to be 
an array of patrol stations fixed in the centers of their 
respective patrol areas. It 1s hypothesized that this two- 
dimensional array 1S an upper bound for the two-dimensional 
emalytic model. 
1. Assumptions 
Assume n points are located within an L by W (L>W) 
area as shown in Figure 6 below. 
Assume that the location of T is uniformly distri- 
buted over the L by W rectangle. The n ships are to maintain 
station at the center of their respective L//n by W//n 


areas. 
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EequmenO.m=eraxeq Lattice Of Points 


2% Development 


Since each of the n ships is centered in an iden- 
tically shaped area the model developed for any one unit 
fet hold true for all units. The maximum value of D is 
the distance from the center to the corner of a unit's 
patrol area as shown in Figure 7. 

a. Conditional Cumulative Distribution Function 


Define event S(i) as the event Te (region 1) 


with P[S(i)] = l/n, and 


Pos (i) (cyte eb <a iS (2) } 
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Bevel eee) x; of Fixed Lattice Model 


4 
2 
Beet bh = 2(a 0%) 7% = aig? - 1/2 
, Z. 
Pength a = ae , length b = C= + = t/a length c = ae 
2/n 2a 

length d = d, the radius of the expanding circle, 

z = 4, Ja = =e 
9 = 2 cos q = 2cOs ae 
arc length s = dé, are length t = dg 


= area of the segment under s but Outside the rectangle 





= =(ds - ga) = ae cost ( iw ) - (a* - Fy h/2 (_W 
2dvn 2Vvn 
— area of segment under t but outside Ene. rectangle 
2 
= (dt = he) = a? cos = (ae = -(a° - 2) 1/2 
2dv/n 2/n 
Figure 7 
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Refer to Figure (7) and consider the case of 


Piomaccr DeCdmetnucle GlE=badius d Such that, 0 < d < W/2vn, 


2 1 : 
ip UD) CISL) = (Te ) Ta = ————_— =. 
n 
Consider the case of W/2/n < d < L/2vn where 
meme Danaing circle has intersected the longer sides of 
the patrol rectangle, 
are 


(4) = SR - (2a: (i 





*D/S (i) Saya) 


As the circle radius continues to increase, consider 
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peeeeconulatitvye DIStrElbution Function 
To uncondition over the event S(i) one needs to 


take the following summation, 





jah 
ee" 
Poid) = & ) Biscay 
ip 
Since the n conditional cdf's are identical. Therefore it 
1s seen that: 
= : goa (0) a el Ss nae 
5 2/n 
nid - 2nA . W L 
Potd) = ar ae BO epee. St 
5 2/n 2/n 
nid - 2nA. - 2nA, L we 7 1/2 
LW Pee Seow Te an! 
2¥n 


(23) 


Sees vALUATION OF THE APPROPRIATENESS OF TWO-DIMENSIONAL 
BPPROXIMATING MODELS 


Expressions (22) and (23) form hypothetical bounds to 
an undeveloped area patrol model. If the hypothesis is 
correct it would seem reasonable that the two-dimensional 
bounding functions should behave in a manner similar to that 
of the one-dimensional bounding functions. That is one 
would not only expect the fixed array model to yield sto- 


chastically larger results than the Poisson field model but 
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one would also expect to see very close agreement between 

the two models for small values of d. Table 2 which tabu- 
lates the evaluation of the two models for parameters n= 25, 
L = 30000 and W = 20000, indicates that for this case the 
two-dimensional models are indeed well behaved. 

In the one-dimensional case it was found that the average 
of the boundary models was itself a very reasonably estimator 
Sieche analytic model. If the eeSS regarding the two- 
dimensional boundary models is correct then a similar averaging 
process could be considered a better approximation of the 
analytic area patrol model than either boundary model is 
when considered separately. The average of the two-dimensional 


boundary models 1s also tabulated in Table 2. 


Model Gvalues aunts Of Jength) 


0 500 1000 1500 2000 2500 3000 3500 3605 


Poisson field ie 2082 se SS) -408 359 692 ioe “OLS 

fixed point sb 053 aie al OASIS 24 i525 5 SA, 998 0 
array 

average 8 pOD2 oor OS fas - 466 .646 .806 OUo 909 


* for parameters n = 25, L = 30,000, W = 20,000. 


Table 2. Evaluation of Two-dimensional 
approximating models 
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V. ANALYTICAL MODEL WITH RADIUS OF INFLUENCE 


The initial analytic model, equation (12), was developed 
under the assumption that an aircraft carrier be considered 
a point on its patrol line segment. This restriction effec- 
tively eliminated from the previous analysis the known ability 
of a carrier's airwing to exercise considerable influence 
over the environment adjacent to the carrier. We now modify 
the model to allow an examination of the contribution to 
ieee) made by an airwing controlling a radius r on either 
Side of the carrier's position. 

An additional benefit of relaxing the zero radius of 
influence assumption may be a quantitative model of how a 
units' probability of being within a specified distance of 
a random target will vary as the radius of influence is 


changed. 


A. ASSUMPTIONS 

All the assumptions made for the initial analytic model 
remain in effect with the exception of assumption II.A. (5) 
which constrained a carrier to being considered a point on the 
patrol line segment. In this modification each carrier will 
be considered to have a radius of influence, r, on either side 
Sepemewcarrier’s position. If the trouble spot T occurs within 
+ r of the carrier's position the distance from the carrier 


to the trouble spot is defined as zero. Distance to the 
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trouble spot, T, will be measured to the edge of a carrier's 
radius of influence instead of to the aren actual 
position as was done in Chapter II. 

There are no restrictions regarding the movement of a 
carrier along its patrol line segment. This means that the 
second carrier would not be restricted from periodically 
een influence Over a portion of the first r miles of the 


third carrier's line segment (see Reet Picgiimes (Si). 


Bee MODEL DEVELOPMENT 
i Cumulative Distribution Function 
In the initial analytical model, 


- Is 
D = min(D,,D,,---,D_) where De (0,7). 


ime the current modification 


, L 
H = min(H,,H,,-.-,H_) where He (0, ~-r) 
and since D; = A Troe a= be 2s. on Ehen 
D-r 1EOuE ra). <= L 
_ —n 
H = 
0 male 0 IB, jae 


iNew £or h > OQ, 
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FIGURE 8. 


the position of the —. carrier 


mine haaditus Of Intiuence On either side 
of a carrier's position 


length of the line to be patrolled 


the number of carriers patrolling 
the line 


location of a randomly occurring 
trouble spot 


the distance measured from the random 
trouble spot to the nearest extremity 
of the it carrier's radius of influence 


Schematic and Symbology for the Radius of 
Influence Model 
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The radius of influence model is 
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which is the initial analytic model [expression (12) ] 


evaluated at D= her. 


ee eeoenal Rate Of Return With Respect to Radius of 
Influence Model 


When the radius of influence model is differentiated 
with respect to r the result is the equivalent to evaluating 


@a@empart of equation (12) at htr: 
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At r = 0 this equation is the true pdf associated 
with expression (12); at values of r > 0 this function gives 
mie wexpected marginal return in terms of P(H<h) for a unit 


change in the radius of influence. 


C. MODEL EVALUATION: 

The following plots were developed to assist in visualizing 
two specific realizations of the cumulative distribution 
function and marginal rate function for the radius of 


influence model. 


L n UN Ol=l ees epko}s Figure 

10 10 ear 9 

10 nO marginal rate a0. 
30000 6 exhin ai 
30000 6 marginal rate 12 

The parameters L = 30,000 miles and n = 6 ships were 


chosen as being representative of the imaginary but roughly 
realistic situation described in the Introduction. 

mae olot in Figure 11 appears to be guite linear for 
all values of r up to and in some cases beyond the r = Q.0 
median value of 1330 miles. It is also of interest to note 
mae, £Or large ratios of =, that for a given h a constant 
Miemeace if © results in a near constant increase in P(H<h). 
This characteristic is not observed in the L/n = 1 case. 

By the time the distance from the target has reached 


3500 miles the n = 6, L = 30000 cdf has exceeded 0.95. In 
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Sebemewobas, LrOrlall values of r > 0 ID 10s SB Sekeh Be 9 5e. 
While this may be a counter-intuitive result it is not 
unreasonable if one remembers that a value of h = 3500 miles 


encompasses a 7000+ 2r mile portion of the total line. 


D. EMPLOYING THE RADIUS OF INFLUENCE MODEL 

We consider some examples which illustrate use of the 
model. The operational commander may desire to know the 
minimum number of units necessary to patrol a 30000 nautical 
mile coast such that the nearest unit has at least a 603% 
chance of being within a 25 day transit of a random spot 
along the coast. Further he may desire to Know if there 
1S any substantial advantage to loading the carriers with 
Type A aircraft (180 nautical mile radius) as opposed to 


loading with Type B (300 nautical mile radius). The units 


steam at 25 kts. 


i 


Serve equation (25) for n at h E500) .6 = 180, b= 730000 


elamee (i <h) >.6. The result is n + 6. Equation (25), 


i 


evaluated at n = 6, L = 30000, h 1500, r = 180 results 


in 


PE Onoll. > 0.6, 


aS required. 


The marginal rate of return function (26), with the 


above conditions gives 
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£., (h) pm 0-000 3), 


which indicates one could expect a 0.0003 increase in 
FY, (A) for each unit increase in xr. The Type B aircraft's 
300 nautical mile radius will increase Fi, (h) by approximately 
moo — 180) (0.0003) = 0.036 over the value of 0.611. (The 
actual increase by solving expression (25) for both R = 180 
aaa R = 300 is 0.03959.) 

As a result an operational planner would inform the 
commander that it would require a minimum of six carriers 
to patrol the line. Additionally, based on an MOE of 
igo d), there is no substantial advantage to be naa by 


loading a particular type aircraft on the six carriers. 
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COF FOR FIXED VALUES 
RADIUS OF INFLUENCE 
L= 10 N= 10 
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VI. EXTENSIONS OF THE ANALYTIC MODELS 


During the course of the development of this thesis 
several potentially interesting extensions were discussed. 
Some of these ideas are presented below as possible areas 


of future investigation. 


A. MULTI-DIMENSIONAL PATROL AREAS 
The assumption limiting a carrier to patrolling along 
a line is unnecessarily restrictive. A natural extension 
would be to expand the patrol line to a patrol area of 
length L and width W such that each carrier would be assigned 


a L/n x W patrol area as illustrated below: 





T, the trouble spot, would be restricted to either occurring 
along the (0,L) line, or perhaps to occurring anywhere within 
the (L,W) area. 

An extension which moves even farther from the patrol 
line model, but which naturally follows the previous exten- 


Sion is to consider a grid imposed on an open ocean area with 
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one carrier assigned each rectangular patrol area. T would 


occur randomly anywhere within the patrol area boundaries. 





Figure 13. n carriers randomly patrolling a 2 dimensional area. 


The radius of influence of each carrier could also be 


incorporated into the two-dimensional models. 


B. GAME THEORY 

The stipulation that the location of the trouble spot 
1s a random variable with a particular probability distribu- 
tion may be realistic in most Situations: involving a naval 
presence mission. Where an adversary relationship can be 
postulated this stochastic assumption as well as the assump- 
tion of independence between the location of T and the 


positions of the carrier units (X,'s) may be untenable. If 
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the adversary 1S capable of some degree of knowledge of the 
Saprler S positions, through radio direction finding, 
shadowing or satellite reconnaissance, he is capable of 
choosing T so as to maximize his return (perhaps by maxi- 
mizing the expected distance to the patrolling force). [In 
this case the problem takes on a game theoretic aspect, and 
results from game theory might be applicable. 

If, through careful analysis of an adversary's intelli- 
gence gathering capabilities, force levels and historic 
political intentions, one is able to determine an a priori 
mesertbution of the location of T, then it is possible to 
use the probability integral transform to transform the a 
aro at Stri bution into a Uniform (0,1) distribution. This 
might enable one to utilize the models developed herein 
and then transform the results back to the original a 


BrLierl space. 


eee CONOMIC COST OF VARYING THE RADIUS OF INFLUENCE 

There are economic costs involved in varying the radius 
of influence of a carrier. Through the use of the radius 
of influence model developed herein and basic cost considera- 
tions, an analysis of the cost tradeoffs involved in achieving 
Emeatcicular improvement in P({D<d) could be developed. This 
development could include economic tradeoffs in terms of 
weetous combinations of: 

1. Increasing a carrier's radius of influence 


2. Increasing the transit speed of the carrier 
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3. Increasing the number of carriers 
4. Increasing intelligence concerning the location 
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